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Abstract 

We suppose: (1) that the ground-state eigenvalue E = F(v) of the Schrodinger 
Hamiltonian H = —A + vf(x) in one dimension is known for all values of the coupling 
v > 0; and (2) that the potential shape can be expressed in the form f(x) = g(x 2 ), 
where g is monotone increasing and convex. The inversion inequality f(x) < /(^r) 
is established, in which the 'kinetic potential' f(s) is related to the energy function 
F(v) by the transformation: {/(s) = F'(v), s = F(v) -vF'(v)}. As an example / is 
approximately reconstructed from the energy function F for the potential f(x) = 
ax 2 + b/(c + x 2 ). 
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1. Introduction 

We suppose that a discrete eigenvalue E = F(v) of the Schrodinger Hamiltonian 



is known for all sufficiently large values of the coupling parameter v > and we try 
to use this data to reconstruct the potential shape /. The usual 'forward' problem 
would be: given the potential (shape) /, find the energy trajectory F; 'geometric 
spectral inversion' is the inverse of this, that is to say F — > /. 

This problem should be distinguished from the 'inverse problem in the coupling 
constant' discussed, for example, by Chadan and Sabatier [1]. In this latter problem, 
the discrete part of the 'input data' is a set {wj of values of the coupling constant 
that all yield the identical energy eigenvalue E. The index i might typically represent 
the number of nodes in the corresponding eigenfunction. In contrast, for the problem 
discussed in the present paper, i is kept fixed and the input data is the graph (F^v), v), 
where the coupling parameter has any value v > v c , and v c is the critical value of v for 
the support of a discrete eigenvalue with i nodes. There are strong indications on the 
basis of studies involving the inversion of the WKB approximation [2] that inversion 
with a fixed i becomes more efficient as i increases (and the problem becomes more 
classical). However, the present paper will be concerned only with inversion from the 
ground-state energy function F (v) = F(v). 

By making suitable assumptions concerning the class of potential shapes, theo- 
retical progress has already been made with this inversion problem [3-5]. In Ref. [4] 
a 'concentration lemma' is proved. If we suppose that Hip = Eip and ||^|| = 1, this 
lemma quantifies the monotone increase in concentration towards x = of the prob- 
ability density ip 2 (x,v) with increasing v. In Ref. [5] this lemma is used to establish 
the uniqueness of the potential shape / corresponding to a given energy function F. 
The class of potentials for which this uniqueness proof applies are those non-constant 
potential shapes / which are symmetric, continuous at x = 0, piecewise analytic, and 
monotone increasing for x > 0. The 'envelope inversion' discussed in Ref. [5] involved 
a class of potentials that could be expressed as a smooth monotone transformation 
f(x) = g(h(xj) of a soluble potential h(x). The approximation obtained was ad hoc in 
the sense that nothing was known a priori concerning the relationship between the 
approximation and the (unknown) exact potential corresponding to the given energy 
function F(v). In the present paper we establish an inversion inequality for a spe- 
cial case of envelope inversion, namely the case in which the 'envelope basis' is the 
harmonic-oscillator shape h(x) = x 2 . Thus we assume that the potential shape f(x) 
has the representation 



H = -A + vf(x) 



(1.1) 



/(*) = 9(x 2 ) 



(1.2) 
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where g is monotone increasing and convex (g" > 0). This is a strong assumption 
but, as we prove in Section (2), it yields a corresponding strong result, that is to say: 

/(*)</(i)> (1-3) 

where f(s) is the 'kinetic potential' corresponding to the potential f(x). The param- 
eter s is equal to the mean kinetic energy < - A> and, in terms of s, the eigenvalue 
F(v) may be represented [6] exactly by the semi-classical expression: 

E = F(v) =min{s +vf(s)}. (1.4) 

The transformations F <-> / are essentially Legendre transformations [8]. This is 
so because we know [4] that F and / have definite and opposite convexity; more 
particularly, we know 

f"(s)F"(v) = ~ <0. (1.5) 

The transformation in the direction needed here F — > / will be given explicitly in Sec- 
tion (2) below where we also prove the inequality (1.3), the main result of this paper. 
In Section (3) we discuss an example for which we compare the upper approximation 
given by (1.3) with the corresponding exact result. 

2. Proof of the inversion inequality 

We suppose that the exact normalized wave function corresponding to the poten- 
tial vf(x) is given by ij){x,v), where the coupling parameter v > 0. Thus (ip,Hip) = F(v). 
We know how this total expectation value is divided between kinetic and potential 
energies for, in more detail, we have 

< - A> = (V>, -Ai>) = F(v) - vF'(v) = s, 
<f> = ( 1 p,fiP) = F'(v) = f(s). 

These equations also define the kinetic potential f(s) parametrically in terms of the 
parameter v > 0. We first use Heisenberg's uncertainty inequality which gives us 

< - Axx 2 > = s<x 2 > > (2.2) 

4 

We now consider 

f(s) = <f(x)> = <g(x 2 )> > g(<x 2 >). (2.3) 

This inequality follows from Jensen's inequality [7] and the fact that g is convex. By 
applying (2.2) in (2.3) we find 
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Finally by letting x = ^ we establish the inversion inequality 

/(*)</(^)- (^) 

□ 

Since the transformation in the direction F -» / is already expressed by (2.1), 
the upper approximation provided by the inversion inequality is now completely de- 
termined. 

3. An example 

We consider the potential shape given by 



f(x) = ax 2 + b/{c + x 2 ), a,b,c>0. (3.1) 

The case a = b = c = 1 is illustrated in Fig.(l) which shows the potential shape 
f(x), in the inset graph, and also the ground-state energy function F(v) generated 
from it. In Fig. (2) the upper approximation A obtained by the inversion inequality 
is shown along with the exact potential shape / itself. The set of corresponding 
'exact' wave functions tp(x,v) are also shown for 3 x 1CT 4 < v < 10. The wave-function 
normalization is arbitrarily taken here to be ip{0,v) = 20, so that the graphs fit on the 
same figure as the potentials. As the coupling v increases, the wave functions become 
monotonically more concentrated near zero, in agreement with the 'concentration 
lemma' mentioned in Section (1). 

3. Conclusion 

Although the assumption behind the inversion inequality is strong, the fact that 
such an inequality exists may be important, especially if it can eventually be gener- 
alized. The expression of this result in terms of kinetic potentials could be avoided 
in principle. However, the representation of the energy functions F(v) in terms of 
f(s) has already yielded some very effective bounds in the forward direction and it 
is natural to explore this same apparatus for the more difficult inversion problem. 
For example, in the forward direction the envelope method [6] may be expressed 
succinctly as 

f(x)=g(h(x)) f(s)*g(h(s)), (3.1) 

where w=> if g is convex and «=< if g is concave. Once one has such an approx- 
imation for f(s), it can immediately be inserted in the expression (1.4) to yield an 
approximation for the corresponding eigenvalue E = F(v). In the present paper we 
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have found one case h(x) = x 2 for which an inequality is retained for the inverse 
problem. In Ref. [5] we also explored the idea of inverting the Rayleigh-Ritz varia- 
tional method and we obtained an inversion approximation with respect to a chosen 
family of 'trial' functions. However, unlike the situation in the forward direction, the 
inversion approximation obtained was again not an inequality. Our experience with 
this problem so far suggests that it is difficult to generate potential inequalities for 
geometric spectral inversion. 
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Figure (1) The potential f(x) = x 2 + 1/(1 + x 2 ) is shown in the inset graph, along 
with the corresponding ground-state energy function E = F(v). The aim of geometric 
spectral inversion is to reconstruct / from F. 
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Figure (2) The approximation A obtained from the inversion inequality (2.5) is 
compared to the exact potential /. The family of corresponding 'exact' wave functions 
tp(x,v) satisfying ip(0,v) = 20 is also shown: the wave functions become monotonically 
concentrated towards zero as v is increased from n = 3x lCT 4 to 10. 



